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MULTIPLE  ELECTROMAGNETIC  SCATTERING  FROM  A  CLUSTER  OF  SPHERES. 

VOLUME  II.  SYMMETRIZATION 

1.  INTRODUCTION 

In  a  preceding  paper  *  (hereafter  referred  to  as  Volume  I),  we 
proposed  a  method  to  study  the  scattering  properties  of  a  cluster  of 
spheres.  Our  approach,  aimed  to  improve  Rayleigh-Debye  theory  for  the 
scattering  of  electromagnetic  waves  from  molecules,  is  quite  general  as  it 
does  not  imply  any  restriction  either  on  the  effective  Index  of  refraction 
or  on  the  size  and  geometry  of  the  cluster.  The  only  approximation 
required  is  the  truncation  of  the  multipolar  expansions  used  to  represent 
both  the  incident  and  the  scattered  waves  as  well  as  the  field  within  the 
spheres.  The  expansion  coefficients  of  the  scattered  wave  are  the  solution 
of  a  system  of  linear  nonhomogeneous  equations  whose  order,  for  large 
clusters  and  for  convergency  reasons,  may  be  rather  high.  If  the  cluster 
possesses  symmetry  properties,  however,  as  is  the  case  for  actual 
molecules,  group  theory  can  be  used  to  get  the  above  system  in  factorized 

Ji 

form.  This  is  easily  done  by  expanding  both  the  Incident  and  the 
scattered  fields  in  terms  of  linear  combinations  of  multipoles  transforming 
according  to  the  rows  of  the  irreducible  representations  of  the  symmetry 
group.  The  symmetr ization  need  not  be  applied  to  the  field  within  the 
spheres  for,  as.  seen  in  section  4  of  Volume  I,  the  expansion  coefficients 
of  the  internal  field  can  be  eliminated  from  the  final  equations. 

Although  the  techniques  we  are  going  to  use  are  well  known  in  molecular 
physics,  two  fundamental  differences  should  be  borne  in  mind.  First,  the 
electric  and  magnetic  multipoles  with  the  same  L  have  opposite  parity  so 
that  they  transform  differently  under  improper  rotations.^  Second,  the 
inhomogeneities  of  the  factorized  systems  depend  on  the  row  index  of  the 
multidimensional  irreducible  representations.  As  a  consequence,  we  have  to 
solve  all  the  systems  arising  from  the  factorization  procedure.  It  will 
become  apparent,  however,  that  these  peculiarities  of  the  parent  problem  do 
not  cancel  the  undoubted  advantages  of  group  theoretical  techniques. 
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2.  SYMMETRY  ADAPTED  MULTIPOLAR  EXPANSIONS  OF  THE  FIELDS 

The  system  whose  scattering  properties  to  be  studied  Is  a  cluster  of  N 
nonmagnetic  spheres  whose  centers  lie  at  Rq  end  whose  red  1 1  and  (possibly 
complex)  refractive  Indexes  are  b  and  n  ,  respectively.  The  cluster  Is 

Cl  ft 

referred  to  a  fixed  system  of  axes  with  its  origin  at  the  center  of 
symnetry,  and  the  direction  of  incidence  of  the  Incoming  plane  wave 
determined  by  the  direction  cosines  of  its  wave  vector. 

Let  us  assume  that  the  cluster  Is  left  unchanged  by  the  operations  of  a 
group  G  of  order  g.  The  effect  of  any  of  the  operations  is  to  permute, 
among  themselves,  the  spheres  In  the  cluster,  but  in  general,  not  all  the 
spheres  are  linked  to  each  other  by  a  group  operation.  Therefore,  we 
partition  the  cluster  In  sets  of  spheres  which  are  transformed  Into  each 
other  by  the  operations  of  G.  He  remark  that  this  partitioning  does  not 
imply  renumbering  the  spheres,  but  only  that  we  associate  to  each  site 
Index,  a  ,  the  appropriate  set  Index,  o  .  The  transformation  properties  of 
the  multipoint  centered  at  the  sites  of  the  o-th  set  are  then  easily  found. 
Indeed,  lf/L(kr)  denotes  a  spherical  Bswsl  or  Hankel  function,  then 


where  the  vector  spherical  harmonics,  X^,,  are  defined  according  to 
Jackson*  and  the  Irreducible  spherical  tensors,  T* ,  according  to 
Rose*'.  The  ebove  equations  show  that  both  the  magnetic,  equation  (la) 
and  the  electric,  equation  (1b),  2L-po1es  transform  under  ^proper" 
rotations  according  to  the  representation  of  the  full  rotation 

group. 7  Accordingly,  let  S  be  an  operation  of  G  such  that 
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with  a  «nd  0  In  the  tame  sat,  of  course,  and  let  0$  be  the  associated 
operator.  If  S  Is  a  "proper"  rotation,  then® 

0f/L(kro)XLM  (r0)  - /L  (k»fe  (rB)  (2a) 


0svx/L(kro)XLM(ra)  vxMkrgjD^SjXu,.  (ft)* 


(2b) 


If  S  Is  an  " improper "  rotation,  however,  the  magnetic  and  electric 
2L-poles  have  opposite  parity,  so  that  the  right-hand  side  of  equation 
(2a)  must  be  multiplied  by  (-)L  end  that  of  equation  (2b)  by 
(-)L+1  ,  while  the  argument  of  Is  to  be  understood  as  the  proper, 

rotational  part  of  S  itself.  Therefore,  if  the  joint  group  G  also  includes 
improper  rotations,  to  obtain  the  symmetrized  combinations  of  multipoint 
belonging  to  the  rows  of  the  v-th  irreducible  representation  of  dimension 
gv  we  have  to  apply  the  projection  operator s7» 9, 10 

'i  •  tEC‘(SI0,  (3) 

’  s 


both  to  the  magnetic  and  the  electric  2L- poles.  Accordingly  we  write 


HVP°  -  £  £  avpah  (kr  )X,  (Pj 
~NL  oTomTn  NLM  L  a  -LM  'a 


(4a) 


and 


Kvpo 

'  N*L 


£  £  *C5f 

a €o  m€n*  ^  * 


vxht(kra)XiJ||(J#) 


(4b) 


for  the  combinations  of  magnetic  and  electric  2L-poies  centered  at  the 
sites  of  the  o-th  set.  The  superscripts  v ,  p  indicate  that  the  combination 
belongs  to  the  p-th  row  of  the  v-th  irreducible  representation,  and  the 
index  N  recalls,  when  appropriate,  that  one  can  get  more  than  a  set  of 
basis  functions  for  a  given  L.  Of  course,  If  G  includes  only  proper 
rotations  -  *nlm  *nd  ^hl  ~  *  The  *cattered  by  the 

cluster,  equation  (3)  of  Volume  I,  can  be  written  In  symmetrized  form  as 


f* 


(!) 


the  corresponding  expression  for  IB*  boing  obtained  through  the  Maxwell 

l 

equation  IB  =  ^pxE;  the  Index  q  ■  il  according  to  the  polarization.  It  is 
useful  to  notice  that  since  we  work  with  unitary  irreducible 
representations,  as  shown  by  the  structure  of  the  projection  operators, 
equation  (3)  ,  the  coefficients  end  hsve  the  property 


(6) 


Now  we  have  to  decompose  the  incident  field,  equation  (l)  of  Volume 
I,  Into  parts  belonging  to  the  rows  of  the  irreducible  representation  of 
G.  This  is  easily  done  for  the  projection  operators,  equation  (3),  has 
the  completeness  property? 


v,P 


PP 


«  1 


(7) 


?: 


Therefore,  we  can  write 

E  Knu|(k)JJ 

E  £  H„u.w[ 

VP  LM  l 


JVP  ♦  nLvp 

LM  t  LM 


<kr)XLM(r)  ♦  n^x.Jc(lcr)XCM(r) J 

1 


(8) 


wher*  “nu. 


is  defined  by  equation  (2)  of  Volume  I  and 

*  M’ 


■§C«»JL<kr)?Ur<C) 

(9a) 

M'  * 

(9b) 

The  analogous  equation  for  1B^  is  easily  obtained  through 

the  relation^ 

._«>  _  <» 

1Bn  »  nEn 

(10) 

Of  course,  even  in  this  case,  if  G  does  not  include  "improper"  rotations 


rvp  -  dvp 

LMM'  LMM* 


-  ■  HZ 
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anna 


This  completes  the  symmetrizetlon  procedure  for,  es  noticed  in  Volume 
I,  the  interne!  field  need  not  be  symmetrized  and  is  therefore  given  by 
equation  (4)  of  Volume  I,  which  Is  shown  below  for  convenience 

§r  •  Z  [cxtr.iur.)  ♦  <"•> 

'C  ■sltXW'CXwy]  (,,b> 

where  R®  and  S®  are  still  the  solutions  of  equation  (5)  of  Volume  I. 


3.  EQUATIONS  FOR  THE  COEFFICIENTS 

The  equations  for  the  coefficients  4vp°  and  Bvp°  of  the  scattered  wave 

T)NL  hnl 

can  now  be  easily  found  through  the  same  procedure  as  used  in  Volume  I. 
First  of  all,  we  rewrite  equations  (5)  and  (8)  in  terms  of  multipoles 
centered  at  the  single  site  ,  by  means  of  the  appropriate  addition 
theorem1 1 

11  =  V  vl  V  >ivP0  r  *VP«I,  tUm  \Y  tl  \  x  V*  Bvpcr  ^Kvpo  1„  .  /.  \ v  /*  X 


E  =  E  E  E^P°  E  »vpah  (kr  )X  (r  )  +  Z  Bvpa  Ebvp°  c*xh  (kr  )X  (r  ) 
v^olItT  nNLM6NNLML  ®-u*-a  7  ^L«TN*N'L"r  L  “  ?Lfil  “ 

+?  ?  4iT£  jL-  (kra^  LM-  (U  +  *  WiM  ^XjL.  ^  LIT  (^0 

+  +ClmHL.  fcr«>!w(^a)]|««) 

"  S  LM'  JL“(kro}'L“M“  (>a)  +  L  L--M-LM-  ^XjL-  H"M" 

^  ^  nd‘^M'  ^  L"M“Ur  jL"^kr0^L“M"  +  *^L"M"LM'  ^r0)-L“M^-0)]j(l3) 


W  =  y*  U 

*1  LM  HLM 


E  E 

vp  M* 


Analogous  equations  hold  for  IB^8’  and  1B^  while  the  matrix  elements  of 
the  dyadic  Greer.'*  function  ^  end  ere 

defined  by  equation  (8)  of  Volume  I.  Moreover,  we  recall  that  equations 
(6)  and  (7)  are  valid  only  for  ra<  Rag  «  Jr^  -  R  I  ;  i.e.,  near  the 


surface  of  the  a-th  sphere.  How  we  proceed  as  in  Volume  I;  l.e. ,  we  take 

A  A 

the  dot  product  of  equations  (11),  (12),  and  (13)  in  turn  *ithra^ra), 

X?  (r  )  and  r  xX*  (r  )  and  get  the  radial  and  tangential  components  of  the 

fields.  Imposition  of  the  boundary  conditions  and  integration  over  the 

angles  then  yields,  for  each  v  ,  p  ,  a  ,  six  equations  from  which  the 

coefficients  of  the  Internal  fields,  C° 

nt»" 


and  D  ,  can  easily  be 

ntm’ 


eliminated.  This  circumstance  clarifies  the  Inessentiality  of  the 
symmetr ization  of  the  internal  fields,  and  yields,  for  each  v  ,  p  ,  a  ,  two 
equations  involving  only  the  4's  and  B1 s  as  unknowns: 

rrrirEf*  «  «  f*,r.*aS  W-i’""*!!:  . bv|>i s''1” 

.  It'JSA  LtMmOxl  Lj  tm  LM  /  NLM  HNL  *£2*  tm  LM  N'LM  HN’L 

T  DcT  L  *  N  MtN 


EE-  (jc)[ncvp  t* 

LM  M'  HLM  *  L  LMM*  £mLM*  LMM 


1 

r  fcmLM'J 

E E E  j  E E  [•! Y'Ca  )  CC  * SjE  * 

T  Bex  L  N'  M€N*  N  M€N 


( 1  Aa) 


EE- 

LM  M* 


HLM 


(s)rcvp 

L  LMM*  tmLM' 


nd1*  1 

LMM'  fcwLM'J 


aB  vpB  vpx 
a  A 

tm  LM  NLM  n  NL 


(14b) 


where  the  quantities  and  5°  are  defined  by  equation  (10)  of  Volume  I. 
Now,  since  the  spheres  within  a  set  are  identical  to  each  other,  the 
quantities  /?£  and  S®  are  actually  independent  of  the  site  index,  a  ,  but 
dependent  on  the  set  index,  o.  This  circumstance,  together  with  the  ortho¬ 
gonality  relations,  equation  (6)  ,  allows  putting  equation  (14)  in  a  more 
symmetrical  form.  Indeed,  multiplication  of  equation  (14a)  by(  avf°)*and 

rptm 

of  equation  (14b)  by  (bv?°)*  and  summation  over  the  a 's  belonging  to  end 
over  m,  yield,  for  each  v,  p,T  ,  the  equations 


E  (•„  kl"  **l  w)- * +E  £  ("••)  ®vp  ■  -pvp 

$  ' n  L  *J  n  /  n*  ~  r»  n»*  nr 


(15a) 


e(«Js.1-‘*bV  (•>Vp+E*v<..< 

r  \ r  *  L  »J  rv  /  n«*  8  r‘* 


.  vp 
)  A  »  -Q 
n» 


vp 

nr’ 


(15b) 


12 


In  equation  (15)  we  put,  for  the  sake  of  simplicity,  r  s  (o.t.n)  , 
s  s  (t,L,N)  r's  (o,  i  ,  n'), s’  a  (t ,  L  ,  N* )  end  define 


#“(->  *EE(avp>*°6  a'jpe 

am  0M  nim  *"«■*  n«-M 

*>,.)«  ££(avP°)*Kof}  bvp6 

"  am  6M  nlm  NLM 


(16) 

(17) 


with  an  obvious  meaning  of  the  parameters  e,m.  The  quantities  ffV((e)  and 
X^(e,m)  are  identical  to  (m)  and  K^(m,e),  respectively,  but  for  the 
mutual  exchange  of  the  a's  with  the  b's.  Moreover 

P'">  £E£„  (k)  (nc'-P  l°  ♦d'-P  j‘  1 

nr  LMonM’  nLM  L  LMM'  £mLM’  LMM'  ftmLM'J  V  ntm  /  *  B#' 


Qvp  W  (fe)  |cvp  La  +  ndvp  «Ja  I  /bvp“V 

n*'  LM  am  M-  nLM  L  LMM'  tmLM’  LMM*  y  trim/ 


(18b) 


We  remark  that  on  the  left-hand  side  of  equations  (16)  and  (17)  the 

superscripts  p  on  fl^and  Kpt  is  missing.  As  will  be  shown  later,  these 

quantities  are  actually  independent  of  the  row  index,  p,  which  has 

accordingly  been  dropped.  Anyway,  the  systems  composed  of  the  equation 

(15)  with  a  given  v,  p  for  all  values  of  o  solve  completely  our  scattering 

problem.  In  fact,  it  was  shown,  in  section  4  of  Volume  I  how  to  compute 

the  relevant  properties  of  the  cluster  from  the  knowledge  of  Aa  and  Ba  . 

*  n  lm  nun' 

the  coefficients  of  the  scattered  field.  Now,  direct  comparison  of 

equation  (5)  with  equation  (3)  of  Volume  I  yields 


a*  ■  £  E  ,ivpVpa  ,  B°  -ZZnupVpa 

HLM  vp  N  nNL  NLM  HLM  vp  N  fgtM 


(19) 


Therefore,  all  the  properties  of  the  cluster  can  also  be  expressed  in  terms 
of  4vp  and  £VP  ,  the  coefficients  of  the  symmetrized  expansion. 

T)NL  T1 NL 

4.  DISCUSSION 

The  consequences  of  the  factorization  effected  by  group  theory  can  be 
best  discussed  by  writing  the  system  with  a  given  v,  p,  in  matrix  form: 


R~l  +  ffv(m) 


#v(e,m) 


Xv(m,e) 


S~l  +  ffv(e) 


AVP 

Pvp 

• 

*VP 

QVP 

(20) 
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Since  equation  (20)  has  the  same  overall  structure  of  equation  (15)  of 
Volume  I,  the  matrix  on  its  left-hand  side  is  easily  recognized  as  the 
inverse  of  the  electromagnetic  T-matrlx  for  the  whole  cluster.12  The 
lack  of  spherical  symmetry  of  the  cluster  does  not  permit  getting  a 
diagonal  T-matrlx,  but  equation  (20)  shows  that  group  theory  effects  the 
decomposition  of  the  whole  scattering  process  into  "modes  of  scattering" 
belonging  to  the  irreducible  representations  of  the  symmetry  group,  G.  Of 
course,  the  dependence  on  the  row  index  of  the  inhomogeneity  of  equation 
(20)  forces  us  to  solve  all  the  systems  arising  from  the  factorization 
procedure,  while  In  the  case  of  secular  determinants  one  needs  to  solve 
only  one  system  for  irreducible  representation.^  However,  as  explicitly 
indicated  by  the  omission  of  the  superscript  p,  the  T-matrix  for  the  v  -th 
mode  of  scattering  does  not  depend  on  p,  as  we  shall  show  presently.  The 
p- i ndependence  of  R  and  S  follows  from  the  very  definition  of  their  matrix 
elements,  equation  (10)  of  Volume  I.  These  diagonal  matrices  occur  even 
in  the  theory  of  scattering  from  a  single  sphere  and  are  the 
electromagnetic  analog  of  the  transition  matrix  of  quantum  scattering 
theory.1*  The  elements  H^,t(m),  ff^(e),  (m,e)  and  J$,  (e,m)  of  the 

matrices  appearing  In  equation  (20)  ,  are  in  turn  the  symmetrized 
counterparts  of  and  which,  as  shown  in  the  appendix  of  Volume  I, 

are  the  matrix  elements,  in  the  site  and  angular  momentum  representation, 
of  G,  the  dyadic  Green's  function  for  free  space  propagation  of  spherical 
vector  waves.  The  p- independence  of  #v(m),  ffv(e),  Xv(m,e)  and  Xv(e,m)  is 
then  a  direct  consequence  of  the  invariance  of  G  under  the  symmetry 
operations.  Therefore  the  whole  T-matrix  for  the  v-th  model  of  scattering 
turns  out  to  be  independent  of  the  row  index,  p.  This  circumstance  greatly 
reduces  the  computational  work  in  the  case  of  multidimensional  irreducible 
representations.  Of  course  the  size  of  the  factorized  systems  depends  not 
only  on  the  number  of  spheres  in  the  cluster  and  on  the  number  of  L-values 
included  in  the  multipolar  expansion  of  the  scattered  wave,  but  also  on  the 
structure  of  the  symmetry  group.  Therefore  the  effect  of  the  factorization 
can  be  Illustrated  only  through  examples.  Table  1  shows  the  order  of  the 
systems  to  be  solved  for  a  cluster  of  5  spheres  with  point  group  T^  the 
CHj^  molecule  and  the  SO^*  ion  have  Just  this  structure).  Although 


tabic  1  considers  tha  case  In  which  tanas  up  to  L  •  4  ara  Included  In  the 
multipolar  expansions,  from  tha  discussion  of  section  6  of  Volume  lf  wa 
expect  that  terms,  up  to  and  Including  L  ■  3  will  be  quite  sufficient  to 
gat  fairly  converged  values  for  tha  scattered  field,  even  If  the  na's  are 
not  close  to  unity. The  usefulness  of  group  theory  requires  no 
further  comment. 


Table  1. 


Lh 

I 

2 

3 

4 

A1 

1 

2 

6 

10 

*2 

1 

2 

6 

10 

E 

2 

8 

12 

20 

Fl 

4 

10 

19 

30 

F2 

4 

10 

19 

30 

U 

30 

80 

150 

240 

Dimension  of  the  symmetrized  and  unsymmetrized  systems  for  up  to  4 
for  a  cluster  of  5  spheres  with  Joint  group  T^.  The  entry  U  means 
"unsymmetrized"  while  the  other  entries  Indicate  Irreducible 

representations.^ 
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